arXiv: 1506.02965vl [hep-th] 9Jun2015 


Avoidance of singularities in asymptotically safe Quantum Einstein Gravity 

Georgios Kofinas^B and Vasilios Zarikas 2,3, 0 
1 Research Group of Geometry, Dynamical Systems and Cosmology, 

Department of Information and Communication Systems Engineering 
University of the Aegean, Karlovassi 83200, Samos, Greece 
2 Department of Electrical Engineering, Theory Division 
ATEI of Central Greece, 35100 Lamia, Greece 
3 Department of Physics, Aristotle University of Thessaloniki, 5fl2f Thessaloniki, Greece 

New general spherically symmetric solutions have been derived with a cosmological “constant” 

A as a source. This A field is not constant but it satisfies the properties of the asymptotically 
safe gravity at the ultraviolet fixed point. The importance of these solutions comes from the fact 
that they describe the near to the centre region of black hole spacetimes as this is modified by 
the Renormalization Group scaling behaviour of the fields. The consistent set of held equations 
which respect the Bianchi identities is derived and solved. One of the solutions (with conventional 
sign of temporal-radial metric components) is timelike geodesically complete, and although there is 
still a curvature divergent origin, this is never approachable by an infalling massive particle which 
is reflected at a finite distance due to the repulsive origin. Another family of solutions (of both 
signatures) range from a finite radius outwards, they cannot be extended to the centre of spherical 
symmetry, and the curvature invariants are finite at the minimum radius. 


I. INTRODUCTION 

The hope that spacetime singularities in classical General Relativity are just an artifact of the high degree of 
symmetries assumed in the solutions was soon lost when the singularity theorems of Hawking and Penrose appeared. 
It is commonly believed that the appearance of spacetime singularities is an indication of the limit of validity of 
General Relativity. A crucial aspect of any successful attempt to quantise gravity is the ability to resolve somehow 
the break down of physics in these singular regions of spacetime. 

Although a considerable effort has been devoted to the discovery of a consistent fundamental theory of nature 
that quantizes gravity, still this quest is one of the most challenging open problems in Science. The perturbative 
quantization of General Relativity leads to a non-renormalizable quantum theory [D-Q. Loop Quantum Gravity 
and spin foam models assume that the spacetime is intrinsically discrete which offer an advance to the problem of 
singularities. Concerning models on continuous spacetimes, one line of thought like String theory suggests to search 
for new symmetries and new degrees of freedom in order to succeed an Ultra Violet (UV) completion of General 
Relativity. Another research path is to keep the same fields and symmetries from General Relativity and work on the 
belief that the associated gravitational action is a fundamental part of the complete theory at the non-perturbative 
level. This is the idea behind the Asymptotic Safety scenario j3-Q| which was first proposed by Weinberg [7;, Q. The 
key issue is the existence of a non-Gaussian fixed point (NGFP) of the renormalization group (RG) flow for gravity. 
If this NGFP determines the behavior of the theory at the UV it can be proved that all measured quantities are free 
from nonphysical divergences. 

The mathematical technique that gave boost to the asymptotic safety scenario is the functional renormalization 
group equation for gravity |9| which enabled the detailed analysis of the gravitational RG flow at the non-perturbative 
level |9l-l34|. This technique uses a Wilsonian RG flow on a space of all theories that consists of all diffeomorphism 
invariant functionals of the metric g^. The model emerging from this construction is called Quantum Einstein Gravity 
(QEG). QEG does not emerge from a direct quantization of classical General Relativity; its bare action corresponds 
to a non-trivial fixed point of the RG flow and is a prediction assuming asymptotic safety. In the asymptotic safety 
program one works with the gravitational effective average action which includes only the effect of the quantum 
fluctuations with momenta p 2 > k 2 , thus r k expresses an approximate description of physics at the momentum scale 
p 2 « k 2 . This effective description helps the development of the phenomenological studies of asymptotic safety. 

The phenomenolo gy o f the asymptotic safety proposal concerning the study of RG improved cosmologies has been 
first appeared in f35U3(| (for a review see IH). The first RG improved solutions of black holes arising from asymptotic 
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safety scenario were derived in [H, (39|. Consequently, other works extended this idea to the Vaidya metric [dO[ , while 
the modified Kerr metric has been presented in |4j (42J . An analysis of the thermodynamic properties of these black 
holes can be found in Phi , while a study of the possibility to include higher derivative terms in the effective average 
action has been presented in Q . New ideas concerning RG improvements have been proposed in lul - lbll . 

There are also attempts to investigate the nature of the microscopic structure of the asymptotically safe quantum 
spacetime [HU-fEU. These works reveal that quantum corrections at high energies (near the nontrivial UV fixed points) 
modify drastically the classical picture. In this direction, but from another perspective, is the study of the possibility 
to have singularity free gravitational collapse [56[. 

The truncated RG flow equations leave two running couplings (with respect to energy), the gravitational constant 
G(k) and the (positive) cosmological constant A (k). Near the non-Gaussian UV fixed point the field G is known to 
approach zero, while on the other hand, the field A scales as the square of the momenta. The scope of our work is to 
investigate the structure of spacetime at high energies and in particular the spherically symmetric configurations in 
the asymptotic safety scenario in the vicinity of the center. This is why we have chosen to study the modified Einstein 
vacuum equations where only A (k) appears. In the same direction is the observation [57, 58] that the cosmological 
constant seems to play a significant role in the understanding of the short distance behaviour. In general, A or G 
energy dependent fields on the right hand side of Einstein equations violate the Bianchi identities. In [59j, certain 
cases with consistent Bianchi identities were presented and solutions were found, however, the case of our interest 
here with A a field unrelated to G was not considered. In several attempts of RG modified black hole solutions there 
is a violation of the Bianchi identities. 

Here, we first find the appropriate covariant kinetic terms that support an arbitrary source field A (fc) without any 
symmetry assumption and construct modified Einstein equations which respect the Bianchi identities. Assuming the 
scaling of A (k) proposed by the asymptotic safety scenario at high energies, we write down and solve exactly the 
consistent quantum gravity improved equations that determine the spherically symmetric solutions. A set of different 
branches of solutions are obtained. The study of their curvature invariants and their geodesic equations reveals a 
variety of smoothening behaviours in clear contrast to the singularities accompanying classical General Relativity. 
Which of these solutions is the physically acceptable one will be determined by the demand to match smoothly with 
the classical Schwarzschild solution at larger length scales. 


II. FIELD EQUATIONS WITH A VARYING COSMOLOGICAL CONSTANT 


In this section we add a spacetime-dependent cosmological constant A(a;) in the 4-dimensional vacuum Einstein 
equations so that G^ = —kg pl ,. Although in the scenario of asymptotic safety, A(fc) is supposed to be determined 
uniquely from the RG flow equations, however, we want at this stage to construct a formulation where A(x) can be 
kept arbitrary (later we will specify the function A(fc)). Obviously, the above equation is inconsistent since it does 
not satisfy the Bianchi identities. We will add an energy-momentum tensor d^y to support A(x), so that the Bianchi 
identities are satisfied. The new equations of motion become 


G 




kg[ii/ T d fj,i /, 


( 2 . 1 ) 


or equivalently 


Rfiu — A g^i, -f- $[iv gfib- ; 


( 2 . 2 ) 


where d = d M „. The Bianchi identities take the form 

H- 

V" = A ;" > (2- 3 ) 

where ; denotes covariant differentiation with respect to the Christoffel connection of g Moreover, obviously d M „ 
must vanish for A constant since in this case the Bianchi identities of m are trivially satisfied and there is no need 
to add anything. Therefore, since the solution of the homogeneous equation d^,,’^ = 0 must be d^ = 0, a partial 
solution of the inhomogeneous equation (12.311 will be the general one. 

It turns out to be very convenient to introduce the field by 

A = A , (2.4) 


where A is an arbitrary constant reference value. We will construct such a tensor d^ from if) and its first and second 
derivatives. We make the ansatz 


d M y = + B(ip)g pv i/j’ p ilj ;p + C + E W + F(ip)g p „ > 


(2.5) 
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where dip = ip.A^- Then, 


i? = (A + 4B)V^W + (C + 4 E)dip + 4F , 


( 2 . 6 ) 


R ij. v = Aip.^ip. v - (77 + B^g^ip’ p ip. p + Cip w „ - + E^jg^dip + (Ae^ - F)g fil/ . (2.7) 

Since d{ip- u ) = (□'0);^ + Rvaip' a , it is from (12.511 

= (A'+B')^;^ + (A+F'pV’Vv + (F' +A+2Fp ; P> ;M;i/ + (C+E)(Op) ;v + FV ; „ + CR^ip* , (2.8) 

where a prime denotes differentiation with respect to ip. Using & equation (12.81) becomes 

pp = (^Fp^AC-Fu)^;/^ + (F'+A-CF-ic 2 )^ W + (C'+A+2B+C 2 )iP' p iP ;ii , v 

+(C+E){\3ip). l/ + (F'—CF+CAe^)ip ;v . (2.9) 

The consistency condition (12.31) becomes 

(^A'+B'+^AC-BC^iP,^ + (e' + A-CE—C^UiPiP^ + (C'+A+2B+C 2 )iP'’ p iP^ v 
+(C+E)(O p). v + [F'-CF+(C - 1)AP] ip- v = 0 . (2.10) 

Since (12.101) must be satisfied for any ip, the various prefactors must vanish separately, giving the following system of 
differential equations 

A' + F' + iAC-FC^O , E'+A-CE-^C 2 = 0 , C" + A+2F + C 2 = 0 

C+E = 0 , F' — CF+(C — 1)AP = 0. (2.11) 


It can be seen that the first equation is redundant and the system gets the form 

A + B = -^C 2 , C' = ic 2 + A , E = -C , F' — CF+(C — l)Ae^ = 0. 


( 2 . 12 ) 


The last eauation in (12.121) for F' is a linear differential eciuation with solution F = c [cr + Af(l-C)e^-f c ], 
where cr is integration constant. Since for ip constant, ■$ AU , must vanish, it arises from (12.51) that for ip constant the 
factor F must vanish. From the form of the previous solution for F, for this to happen for any constant ip, it needs 
to be C = 1 and furthermore a = 0. Thus, although (12.121) is a system of four equations for five unknowns, the above 
requirement gives a unique solution. The system (12.121) has therefore the solution 


A = -- 


F = —- , C = 1 , F = —1 , F = 0 . 


Finally tF,, takes the form 


®nv = ~^P-,pip-, v - + V’w*' - 9p.udip . 


(2.13) 


(2.14) 


Equation (12.141) is the general solution of (12.31) and could also be expressed in terms of A. The self-consistent modified 
Einstein equations HU become 


GW — AcA g pl) ^.9/xt'V 7 ' P<i P\p At ip-p-u 9pA-\ip, 


(2.15) 


containing an external field ip(x) which is given by a favored theory. Note that ip does not have its own equation of 
motion. 

If an extra matter energy-momentum tensor T pl , is included (which will not be the case in the present work) the 
field equations become 

Gpv — Ae^g/jv ^ip;pip-,i/ ~^9pv' l P’ p 'ip;p A ip-,p-v 


g^Uip + 8 nGT^ , 


( 2 . 16 ) 
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where the gravitational constant G is in general also spacetime dependent. Since the various ^-kinetic terms along 
with the A field in (12.1611 are by construction covariantly conserved, the Bianchi identities provide that the energy- 
momentum is not separately conserved but only in combination with G, thus 

{GT Mv y v = 0. (2.17) 

If one insisted on the exact energy-momentum conservation, then extra covariant terms constructed out of G would 
be necessary to be added in (12.161) . However, we do not think that it is pertinent to include such an extra assumption. 
Therefore, since the RG flow equations aim to provide energy dependent A (k),G(k), the system of equations (12.161) . 
(12.171) is a consistent system of equations valid at any energy scale in the presence of matter. The consistency closes 
with the connection of the momentum scale k with some length scale and the only limitation is the restricted knowledge 
possessed about the form of the functions A(fc), G(k) as a result of the RG flow equations. 

There is another approach where the quantum corrections of G and A enter not through the improved field equations 

(12.161) . but through an improved action whose variation gives different field equations. However, in this approach the 
full A and G-dependent terms which make the equations consistent have not been found in general. As for equations 

(12.161) it is not evident if they arise from an action or not. 


III. SPHERICAL SYMMETRY 


We are looking for static spherically symmetric solutions. The general metric is of the form 

ds 2 = —A(r)dt 2 + B(r)dr 2 + r 2 (d0 2 + sin 2 0d<j) 2 ). 


(3.1) 


The external field A(x) carries the same symmetries, so it is A = A(r). For the metric (13.11) the non-vanishing 
components of the Einstein tensor G(( are 


G = 


1 

r^B 


1 


G r _ 


c® _ r 

— c J. 


r 2 B 

A' 


1 

-o T 


B' 

Vb 2 

A' 


2 rAB 
B' A' 2 


MB' 


A' 1 


(3.2) 

(3.3) 

(3.4) 


2 rAB 2rB 2 4A 2 B 4AB 2 2 AB ' 

A prime denotes here and in the following differentiation with respect to r. It is 'ip’ p ip- tP = ■%-, Oip = ^1 _|_ AjL _ _|_ 

and the non-vanishing components of ip-^-v are ip-f,t = —4tr — TTg - ) ' t P;0;9 = = rsin B 8 ^ . 


_ 2, 

Therefore, the three independent components tt. rr, 99 of (12.151) 

1- B B' 


are 


A' 

2 rA 


B' 

2rB 


A' 2 

4A 2 


rB 
A!_ 
rA 
A" 

4 AB + 2A 


1 - B 
^2 

MB 


= -A e*B-^-W + W_ r 
4 r 2B * 


= —Ae^B — 


= -A e+B-¥--*- 


3V> /2 2 if>' A'tjj' 

~4 r 2 A~ 

A'tp 1 


2 A 


B'ip' 
2 B 


-r 


(3.5) 

(3.6) 

(3.7) 


The system of equations (13.5I) - (I3.7I) is satisfied by construction for any ip(r). It can be seen that the equation (13.71) is 
redundant, since it arises from (13.51) . (13.61) and this reflects the r—reparametrization invariance of the system (12.151) . 
Indeed, differentiating (13.61) with respect to r, we obtain together with (13.51) an algebraic system for . We solve 

this system for A",i/j" and we also solve (13.61) for A'. Plugging these A”, ip", A' in (13.71) we find an identification. 

At the UV the RG flow equations provide A = A*fc 2 , where A* is a dimensionless constant. Note that this scaling is 
also consistent with dimensional analysis without the introduction of a new energy scale. It is natural to set a relation 
between the energy and length scale of the form |35| 


k = 


i 

D 1 


(3.8) 


where £ is a dimensionless parameter and D = D(r) > 0 is the proper length of a radial curve defined by dt = 
d9 = d<p = 0 connecting a point with coordinate r to the minimum radius of definition r = r m i n (possibly the origin 
Gnin = 0). Thus, 


D(r)= f V\B(pj\ dp 

J Gmir, 


(3.9) 
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and D'(r) = \/\B(r)\. It turns out that ip' = -2%-, %■ = 2%r- Equation (13.51) is written as 


D" i / r‘ \ „ /9 LT 1 
2 W + r( € ~ a Dl) D + 4_ D~r = °’ 


D ’ B 
1 ' - 2 


£>' 1 


(3.10) 


where a = 3 + eA*£ 2 and e is the sign of B(r). Normally it is e = +1, but since the signature of the metric is not 
known close to the centres which are the regions of high energy, we leave e = — 1 as a possibility (for example, the 
interior Scliwarzschild black hole in the same spherical coordinates has e = — 1). As for A however, we will require 
in the following that A and B have the same signs, so there is always a single time dimension. Of course, r m ; n does 
not enter (13.101) . or to say it different, an arbitrary constant could have been added on the right hand side of (13.91) 
without changing (13.101) . This means that in the end, the normalization in order to find D will come by demanding 
D = 0 for r = r mi n . 

Equation (13.101) has the symmetry D —> A D, r -+ Ar, so under this transformation, (13.101) remains invariant. Setting 


D 

u = — 
r 


equation (13.101) becomes 


Under the redefinition 


equation (13.121) becomes 


ru"+2u' 1 

-/-1- 

ru +u r 


+ -fe— (ru'+u) 2 +4 
r \ u J 


ru'+u 1 
ru r 


n u + u ( a \ . , 2 

2-- hie - 9 )(u + it) 2 + 4 

u+u \ u 2 J y ' 


u + u 

where a dot denotes differentiation with respect to r. Defining 

v = u , 

equation (13.141) becomes 


u + u 


■u 


-1 = 0 , 


2v d 
v + u du 


(v + u) + ^e - (v + u) 2 + 4 


v + u 


-1 = 0 . 


Setting 


y = v + u , 

equation (13.161) reduces to the first order equation 

. dY ( a \ ., 4 0 

2 (y- u — + e - -=■ y 3 + -y 2 - y = o. 

du V u z J u 


Under the replacement 


equation (13.301) becomes 


V = 


Y 


= — 2_A + — 


uJdu 2u^ 


u 2 ^ 2u V u 2 


(3.11) 

(3.12) 

(3.13) 

(3.14) 

(3.15) 

(3.16) 

(3.17) 

(3.18) 

(3.19) 

(3.20) 


Equation (13.201) can be converted to an Abel equation and such an equation is not in general solvable. However, 
very luckily, here the coefficients of (13.201) are of a special form and equation (13.201) can be converted to a separable 
equation. Namely, setting 


z = 


\fu 


(»-;) 


(3.21) 
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equation (13.201) becomes 


dz a — 3 


du 2 u 4 2 u 2 


z = ±\/4r + - + c , $ = 1 — v > 


with general solution 


where c is integration constant. Then, 

1 ± v/cw 3 +ert 2 +/3 

2/ =-- • 

u 

It is obvious that it has to be cu 3 +eu 2 +/3 > 0. 

From the above definitions of the various variables, the metric component B(u) = B(r(u)) is found to be 


(3.22) 


(3.23) 


(3.24) 


B(u) = 


eu 


(l± \/cu 3 + eu 2 +0) 2 

In order to find the metric component A , it is convenient to introduce the variable \ instead of A defined by 

A - —p 2 * 

A- B e ■ 

Plugging this expression in equation (13.61) and subtracting from (13.51) we get 

9\ 1 

/'il' " L 


+-)x'=r-^' 2 - 

r J 2 


(3.25) 


(3.26) 


(3.27) 


Since u' = ±(± - u), ip' = it is ip" = + (| - u)^]. Then, equation (13.271) due to (13.201) becomes 

2(1 ~uy) 2 j^ = y+ — — — - . 

au \u z / y u 

Replacing the solution (13.241) in (13.281) we obtain the equation 


2(»W+l-f)£ 

V 3 / du 

Integrating this equation we find 


dx _ 3cu 2 +2eu (cu 3 + ^-) ^cu 3 +eu 2 +1 — 1 


3/ du cu 3 +eu 2 —j 


u(cu 3 + eu 2 — §) 


2 x(u) = Ini? + In .,' " „ , 3 ± J(u). 


\cu 3 +eu 2 -f\ 


where > 0 is integration constant (not to be confused with of section ED and 


3,2 a 


J{u) = 


i(cu 3 +eu 2 -j)^cu 3 +eu 2 + l — - 


zdu. 


(3.28) 


(3.29) 


(3.30) 


(3.31) 


Then, A is found as 


A(u) = ^ (l± i/ctt 3 + eM 2 + l — — 


cu 3 +eu 2 + l — % 


\cu 3 +eu 2 -f 


3. e ±3(u) 


where A(u) = A(r(u)). Finally, the relation between u and r is found from equation rv! = 1 Mv to be 


(3.32) 


ry—du 
r = roe J 1 -“" . 


(3.33) 
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where ro > 0 is integration constant. Using the solution (|3.24l) we find 

- 1 ° P Tl(n) 


{ \ ' u 

r(u) = — e 
u 


(3.34) 


where 


l(n) = [ - 
J i 


du 


i^/cu 3 +eu' 2 + 1-f ’ 

To summarize, the metric in the non-conventional gauge determined by u is given by 


(3.35) 


ds 2 = —A{u)dt 2 + r(uY 


-cu 3 + eu 2 +1 —§ 


du 2 + {d9 2 +sm. 2 6 dcj) 2 ) 


(3.36) 


It is more convenient to write the metric (13.361) in the form 


ds 2 = lZ(u) 2 — A(u)dt 2 + 


ern 


cu 3 +eu 2 +1 —§ 


du 2 + rg(d9 2 + sin 2 dd<ji> 2 ) 


(3.37) 


where 


K(u) = - 
u 


(3.38) 


A(u)=ed(l±Jcu 3 +eu 2 + l-^) 


\\ 3 / \cu 3 +eu 2 — %\ 


(3.39) 


J(v) = J 


3cu 2 + 2eu 


(cu 3 +eu 2 — ^)^cu 3 +eu 2 +1 —f 


du. 


(3.40) 


The relation of r(u) and TZ[u) is r{u) = r 0 lZ(u), while the relation of A(u) and A{u) is A(u) 
It is useful to define the variable 




v = „+-. 

Since u > 0, it should be v > ^. The metric (13.371) is expressed in terms of v as 


(3.41) 


ds 2 = R(v) 2 — A (v)dt 2 + 


c(v 3 +pv+q) 


dv 2 + rp(dd 2 + sin 2 dd</> 2 ) 


(3.42) 


where 


R(v) = 




(3.43) 


A(v) = ed sgn(c) 1± \/c(v 3 +pv+q) 


2 v 3 


A+pv+q e± J(v) 

I v 3 —(— pv ~\~q “I 


(3.44) 


The integrals I(v), J (v) are 


W1 


dv 




(3.45) 


J(v) = / 


3v 2 +p 


dv. 


(v 3 +pv+g— i) y 7 P(v) 


(3.46) 
































where 


P(v) = c(v 3 +pv+q) 


P = ~ 


3c 2 


_ 1 / 2e a\ 

“ c V27 c2 + “3/ ' 


(3.47) 

(3.48) 


Now, it is r(v) = roR(v) and what remains in the next is the calculation of the integrals /(v), J(v) encountered in 
R(v), A(v). From (13.111) . (|3.34l) and since X(zt) = /(v) we get 


D = r 0 [e* 1 ^ - e^] 


where v corresponds to r m i n . 

Note that A(v) from (13.441) is written as 


while from (13.461) it arises that 


A(v) = etf 


dJ 

dv 


(/FV)±1) 2 

l - P ( v )— 1 — 1 | ' 


>±J(v) 


dP 


VP(P-l) dv ' 

Differentiating (13.501) with respect to v and using (13.511) we obtain the simple differential equation 


The solution of (13.521) is 


1 dA _ 1 dP 
A dv P dv 


A(v) = e(P(v) = e(c(v 3 +pv+q ), 


(3.49) 

(3.50) 

(3.51) 

(3.52) 

(3.53) 


where £ > 0 is integration constant. Therefore, although the integral J(v) in (13.461) looks rather complicated and can 
be converted to typical elliptic integrals, however, the specific combination 3v 2 +p = in the numerator of (13.461) 
implies that e ±J ( v ) is expressed as a simple algebraic function of v. Finally, the spacetime metric (13.421) takes the 
form 


ds 2 = R(v) 2 - eC,c(v 3 +pv+q)dt 2 + 


c(v 3 +pv+q) 


dv 2 +r 2 0 (dd 2 +sm 2 dd(l) 2 ) 


(3.54) 


So, the problem has been reduced to the calculation of the integral /(v) determining R(v). The solution (13.541) contains 
three integration constants ro,c, £ but ( can be absorbed in a rescaling of the time t. Using (13.431) . (13.451) . the Ricci 
scalar Ric of the metric (13.54[) can be found to be 


Ric = _^ e ±2/(v) 


6ec 


2 ^ 
rn 


±2/(v) 


± 27 p ( 3v3 + ^ 2 + 5PV+6<? + I ) +2q+ ^ + ~d 
± ^& 2 + 2 pv+ 3 q+e £) + 2 q+ i^ + l ± Tc^. 


(3.55) 

(3.56) 


IV. SOLUTIONS 
A. General solutions 


Up to this point, we have been quite general. Now, we will restrict our interest to the case where the polynomial 
P(v) has three real roots po , Pi,P 2 , he. p 3 +ppe+q = 0, t = 0,1, 2. Then P(v) is decomposed into three simple factors 
without the need to introduce complex roots. The three real roots exist when 4 p 3 + 27 q 2 < 0 (i.e. |< 7 C 3 | < 2/27), or 
equivalently ||f + 9(3 — a)| < which means 0 < e(a — 3) < ^, and therefore 0 < A*£ 2 < ^■ The roots pe are 
given by 


Pi = 


3|c| 



r 1 

27 l (3 \ 

2n£-\ 

COS 

- arccos | 

- o 

r t |c| ’) 



(4.1) 
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It can be seen that — gj^j- < pi < — < P 2 < 55 ^ < po < The variable v is defined in some interval (finite 

or infinite) such that P(v) > 0 (formally the equality should not be included as P(v) is denominator in (13.541) . but 
it more convenient to leave this notation in the following analysis). Depending on the domain of definition of v, 
the integral J(v) may possess an integrable pole v = pe and/or the non-integrable (logarithmic) pole v = ^. There 
are only four cases for the domain where v is defined: (i) c > 0 , v > po, (ii) c > 0 , e = — 1 , —^ < v < P 2 , (iii) 
c < 0, e = 1, p 2 < v < po , (iv) c < 0, e = —1, — A < v < po- In the following analysis we make use of the identities 
P 0 +P 1 +P 2 = 0, P 0 P 1 +P 0 P 2 +P 1 P 2 =p, P 0 P 1 P 2 = —q- 

• Case (i): c > 0, v > po- 
We use the transformation 


x = 


v—Po 

V p2 


(4.2) 


where the lower limit of the integral /(v) in (13.451) is v = po and the upper v. We find (correcting the formula 7, p. 
262 of [13 and the formula 5, p. 58 of [6l| l 


2 

p( 

o 

1 

> 

P2~Pl\ 

P0-P2^f 

~<T 

1 

o 

P^~Tc 

P 2 -P 1 M 

(P2-^)\/c(po-Pl) 

V \\ 

\j V p2 

5 

1 

o 


\j V-p 2 

P°-Tc 

O 

1 


(4.3) 


where 


F(z, to) = 

is the incomplete elliptic integral of the first kind, and 

11 ( 2 , n, to) = 


dx 


0 \J (1 —x 2 )(l —too; 2 ) 


dx 


0 (1 — nx 2 )\J(l — x 2 )(l — mx 2 ) 


(4.4) 


(4.5) 


is the incomplete elliptic integral of the third kind 1 . The above expressions hold for Po 7 ^ 53 , P 2 7 ^ ^ which mean 
c 3 |g| ^ Tf- I n this case, the function /(v) ranges from zero value at the minimum value of v = po up to a finite 
positive value 1^ for v —> + 00 . 

For the upper branch of the various ± signs, the radial function R(v) extends from the origin r = 0 (for v —> v = + 00 ) 
to a finite radius (for v = po). At this finite radius, the coordinate r max is from (13.431) r max = , while the proper 

distance D is from (13.49[) D max = ?'o(l — e _/ °°), normalizing this way the integration constant ro. The Ricci scalar, 
as it is seen from (13.551) , diverges at the origin signaling a curvature singularity there. Additionally, the metric (13.541) 
close to the origin r ~ 0 takes the following form 

dsl — —e(ce~ 3I °° —dt 2 + -e 7 ” —dr 2 + r 2 {dd 2 + sin 2 9 df> 2 ) , (4.6) 

r c ro 

from where it is also seen that the Ricci scalar diverges as R ~ r~ 2 , while R^ V k\R , ivk ‘ X ~ r ~ 6 (for the Schwarzschild 
metric it is R = 0, R i i VK ,\R >lVK ' > ' ~ r -6 ). Close to the origin r ~ 0 the proper distance is D ~ —L=e^ /oo r 3 / 2 , the 
same scaling as Schwarzschild. For e = 1, although the solution is divergent at the origin, we will show in the next 
section that this divergency is not felt by any infalling timelike particle of arbitrary initial velocity. The reason is that 
the particle, due to a seemingly repulsive force close to the origin, stops at some finite distance and then it moves 
outwards. 

For the lower branch of the ± sign, the consistent solution (with D > 0) ranges from a minimum non-zero 
radius r m i n = ■ T ° _e_ for v = v = p 0 to a maximum radius r max = - r ° eR eo ^ for v = go > po- It is 
go = cos{^ arccos [^c 3 (^ — q)]} for e = 1, 0 < A *£ 2 < ^ - 3 or for e = —1, 3 < A *£ 2 < while 
it is go = —^sgn(cg-l)cosh{|arccosh[^c 3 (^ — q)]} for e = 1, 0 < ^ ~ 3 < A *£ 2 < ^ or for e = —1, 


1 In the literature, F(z,m) is denoted by i^arcsinz, y/rn) in [6cL [6lfl and by i^arcsin z\m) in [62t and Wolfram’s Mathematica, while 
fl(z,n,m) is denoted by II(arcsin z,n,y/rn) in [60| , |[6lj| , by II(n; arcsinz\o:) (with m = sin 2 a) in [62ll and by II(n; arcsin z\m) in 
Wolfram’s Mathematica. 
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0 < A *£ 2 < 3 < The maximum proper distance is D max = r 0 [e 7 ^°) — 1]. However, still here, the Ricci scalar 
diverges at the minimum radius, as it is seen from (13.561) . and we do not analyze this solution further. 

• Case (ii): c > 0, e = —1, — ^ < v < p 2 - 
We use the transformation 


x = 


Po~Pi P2~v 
P2 Pi P0 ~ v 


(4.7) 


where the lower limit of the integral /(v) in (13.451) is v and the upper v = p 2 - We find (in agreement with formula 4, 
p. 262 of (60| and formula 6 , p. 65 of 61) 


I(v) 


(po+^)V c (po-pi ) L M/ P‘2 pi Po-y 


F 


(/ 


po—pi p2 y 


P2-Pl \ 

Po~P\) 


P2~P0 / PQ-Pl P2 y P2~Pl P0 + ^ P2~P1 \' 
P2 + ^ W P2~Pl P0-y’ PO-Pl P2 + ^’ PO-Pl'- ’ 


( 4 - 8 ) 

In this case, the function /(v) ranges from +oo at the minimum value of v = — ^ up to the zero value for v = p 2 - 
For the upper branch of the various ± signs and 3 < A *£ 2 < the consistent solution (with D > 0) ranges from a 
minimum non-zero radius 


X min 


r ° r ~i(e 2 ) 

B 2 + 57 


(4.9) 


for v = v = Q 2 = ^ cos{i arccos [^yc 3 ( 4 — q )] + -^} < p 2 to a finite maximum radius r max = for v = p 2 . The 

maximum proper distance D is from (13.491) D max = ro[l — e^ 1 ^ 2 ^]. At the maximum distance, the Ricci scalar is seen 
from (13.561) that diverges. However, this is not important since the solution found is physically realistic close to the 
minimum radius where the cosmological constant scales as k 2 ; at larger distances a deviation from the k 2 law occurs 
and the solution has to be modified there, while at even larger distances the solution should reduce to the low energy 
Schwarzschild (or Schwarzscliild-De Sitter) solution. The important information carried by this solution is the fact 
that the present branch has at the minimum radius finite the Ricci scalar as well as the other curvature invariants. 
At radii close to the minimum, this solution can also be considered as physically realistic and non-singular since the 
curvature invariants are finite. This is a success of the quantum gravity scenario under discussion and the avoidance 
of singularities is a silent hope of quantizing gravity. 

Yet for the upper branch, there is another solution for 0 < A*£ 2 < 3 < gij which extends from the origin r = 0 for 
v = v = — -A to a finite radius r max = for y = p 2 - However, the Ricci scalar diverges at the origin. 

For the lower branch of the ± sign, the solution ranges from a minimum non-zero radius for v = p 2 to a maximum 
radius for v = — The Ricci scalar diverges at the minimum radius. 


• Case (iii): c < 0, e = 1, p 2 < v < po- 
We use the transformation 


x = 


Po~Pi y~p2 
Po~p 2 v— pi 


(4.10) 


where the lower limit of the integral /(v) in (13.451) is v = p 2 and the upper v. We find (in agreement with formula 5, 
p. 262 of [(3(3 and formula 6 , p. 63 of [ 6 l|) 


i(y) 


2 

(Pi--kW c (Pi~Po) 


F 


( / Po Pi y p 2 
vy p 0 -p 2 v pi ’ 


P0~P2 \ 

po-pi' 


P2-pl jW / Po — Pi y~P2 P 1 --& PQ-P2 P0-P2 Y 
P2-53 W P0-P2 y—pi ’ P2 — ^ po-pi po-pi'- 


( 4 -H) 

The above expressions hold for pi,p 2 7 ^ which mean c 3 q ^ Ay. In this case, the function /(v) ranges from zero 
value at the minimum value of v = p 2 up to a finite positive value for v = p 0 . 

For the upper branch of the ± sign, the solution ranges from a minimum non-zero radius for v = po to a maximum 
radius for v = p 2 - The Ricci scalar diverges at the minimum radius. 

For the lower branch of the various ± signs and 0 < A *£ 2 < ^ — 3, the consistent solution (with D > 0) ranges from 
a minimum non-zero radius 


r m in 



e Aeo) 


(4.12) 
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for v = v = po = — £ cos {^ arccos c 3 (q — A)]} < po to a finite maximum radius r max = pJ-A. for v = p 0 . 

The maximum proper distance D is from (13.491) D max = r^e 1 ^ 0 ^ — e I ^ eo ' > ]. This branch is the second branch we meet 
in the analysis of solutions where the curvature invariants at the minimum radius are all finite and note that the 
signature of this metric is the “correct” one. 

Yet for the lower branch, there are another two solutions. The first one for 0 < A*£ 2 < — 3 extends from a minimum 

non-zero radius for v = p 2 to a finite maximum radius for v = Q 2 = — ^ cos{;| arccos [■=^-c 3 (q — A)] + 4^} > P 2 and 
has divergent Ricci scalar at the minimum radius. The second one for 0 < — 3 < A*£ 2 < g|r or 0 < A*£ 2 < ^ < 3 

also extends from a minimum to a maximum radius for v = P 2 and v = po respectively and also has divergent Ricci 
scalar at the minimum radius. 

• Case (iv): c < 0, e = —1, — < v < p 0 - 

We use the transformation 


x = 


Po-y 

P0~P2 


(4.13) 


where the lower limit of the integral /(v) in (13.451) is v and the upper v = po- We find (in agreement with formula 6 , 
p. 262 of 0| and formula 5, p. 61 of 61]) 


J(v) = 


(po + ^)Vc(pi-po) 


: n 


(/ 


PO —V Po~P2 P0~p2 


P0~P2 P 0 + 3 J PO — Pl 


(4.14) 


This case is very similar to the case (ii). The function /(v) ranges from +00 at the minimum value of v = — up to 
the zero value for v = po- 

For the upper branch of the various ± signs and 3 < A *£ 2 < the consistent solution (with D > 0) ranges from a 
minimum non-zero radius 


r-min 


r 0 

So + ^ 


e -i(eo) 


(4.15) 


for v = v = g 0 = —ff- cos {A arccos \^-c 3 (q — A)]} < po to a finite maximum radius ?' max = —for v = p 0 . The 

maximum proper distance D is from (13.491) H max = ro[l — e At the maximum distance, the Ricci scalar is seen 
from (13.561) that diverges. However, at the minimum radius, the Ricci scalar and the other curvature invariants are 
finite and this is the third case of a non-singular solution which does not extend to the origin. 

Yet for the upper branch, there is another solution for 0 < A *£ 2 < 3 < which extends from the origin r = 0 for 
v = v = — to a finite radius r max = —for v = po- However, the Ricci scalar diverges at the origin. 

PO+3^ 

For the lower branch of the ± sign, the solution ranges from a minimum non-zero radius for v = po to a maximum 
radius for v = — J-. The Ricci scalar diverges at the minimum radius. 


B. A special solution 

In this subsection we study a simpler special solution which arises when the integration constant c = 0 and is not 
included in the general solutions derived above. From equation (13.371) we take the solution 


ds 2 = 


lZ(u) 2 — ei9(eu 2 + l — H-—— du 2 + r^dd 2 + sin 2 d dcf 2 ) 

' O' au JL tt 


(4.16) 


where eu 2 + 1 — § > 0. The function 7 Z(u) = is given from (13.381) for e = 1 by 


H(u) = -e^ a “coos- 
u 


u> a , a = A /- - 1 , 


a 


(4.17) 


while for e = — 1 by 


KW = -( 


( \/a 2 — u 2 + a^-k 

, . A 

— 

, u < a , a = \ 

V u J 

V 


a 

3 


(4.18) 
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The Ricci scalar is found to be 


Ric = —~e 
»o 


±2X(u) 


'9 — 2 a , 2it 2 —3d 2 


3e 


± 




(4.19) 


where 1 (u) = A arccos — for e = 1 , while I(u) = A In ... = = —- for e = — 1 . The upper branch for e = 1 and the 

lower branch for e = — 1, 0 < A *£ 2 < 3 involve the centre r = 0 and the curvature invariants diverge at r = 0. The 
lower branch for e = 1 , a < u < y/a 2 + 1 and the upper brance for e = — 1 extend from a non-zero minimum radius 
outwards and the curvature invariants diverge at the minimum radius. The only significant branch is the lower one 
for e = —1, A*f 2 > 3 which is free of divergencies in the curvature invariants. Namely, this solution ranges from a 
minimum non-zero radius for u = y/a 2 — 1 to a finite maximum radius for u = a and at the minimum radius it is 
seen from (]4.19(1 that the Ricci scalar is finite (the same is true for the other curvature invariants). 


V. THE GEODESIC EQUATION 

Let us first write the metric (13.541) in the form 


ds 2 = - 


eCP(v) 


z* 2I ^dt 2 


-( v -fe) 2 


( y -^) 2 ' ' [1±V?W] : 


■ dr 2 + r 2 (d0 2 + sin 2 9 dcjr 2 ). 


(5.1) 


where v is supposed to be implicitly solved in terms of r from (13.431) . This metric arises because the relation of dv 
and dr is 


(dv \ 2 P(v)e ±2/ ( V ) / e \ 4 

(») = rgl LW ( T -s;) • (5 ' 2) 

To find the timelike or null geodesics in the above spacetime, we orient as usually the coordinate system so that the 
geodesic lies in the equatorial plane 9 = j. Due to the presence of the Killing vectors the energy E and 

the angular momentum J of the particle are conserved, E = —gttjjji J = where A is an affine parameter. 

The magnitude of the 4-vector of energy-momentum is given by the rest mass /i of the particle (/r = 0 for massless 
particles), so 17/^yfx" + T 2 = 0 an d therefore 


/ dr \ 2 / r2. 

\9tt9rr\ (^) =£ 2 -eM(/x 2 +^). (5.3) 

For massive particles we get 

\9tt9rr\(^) =E 2 -e\g tt \(l +—), (5.4) 

where r = A/i is the proper time, E = ^ the energy per unit rest mass and J = A the angular momentum per unit 
rest mass. We consider for simplicity radial orbits, so J = 0 and we obtain 



2 


[1 ±/PW] 2 


,fj2 e ±2I(v) 

- C P(y) 


(y-fj 2 


(5.5) 


For e = — 1, equation (15.51) does not set any additional constraint on the allowed range of v (and accordingly of r) 
that the particle can acquire. However, for e = 1, the right hand side of (15.51) should be non-negative, that means 


(P(v)e* 2I M 

(y-Tc ) 2 


< E 2 . 


(5.6) 


As we explained in the previous section, for the upper branch of case (i) with e = 1 the solution extends from the 
origin r = 0 for v —> +oo to a finite radius for v = po- For this solution, the inequality (15.61) means that for any 
E the allowed for the massive particle range of v is from a minimum value of v up to a finite value v = vj which 
saturates the inequality. This means that for the particle, v cannot approach infinity, and thus r cannot approach 
the origin. The minimum distance Df from the origin that the particle reaches is Df = ro[e _/ ^ v/ -*—Although 
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the curvature invariants diverge at the origin, however the origin is inaccessible by any timelike particle. We can also 
find the proper time r tot needed by an infalling timelike particle starting from an initial parameter v = v, in order to 
reach the minimum distance 


Ttot = r 0 



e I Wdv 

( v -^V?( v -^ 2e2/(v) - p ( v ) 


(5.7) 


By definition, the upper value Vf of the integration vanishes the square root in the denominator of the integral 
ED However, this pole is seen to be integrable, therefore the time r to t needed from the particle to reach the 
minimum distance is finite. This means that the particle bounces at the minimum distance and afterwards it moves 
outwards. Thus, the energy-dependent cosmological constant (along with the other terms needed to make the equations 
consistent) acts repulsively and is able to inhibit an infalling mass, while a simple cosmological constant is unable 
to deter the collapse of the particle. 2 Since the particle is reflected, its geodesic is extendible (at least in a region 
of its minimum approach), so the space is what is called t-complete, i.e. it is geodesically complete for the timelike 
geodesics. 

In this spacetime however, a radially infalling massless particle reaches the origin r = 0 in a finite affine parameter. 
Indeed, we have 



E 2 e 2/ ( v ) 


[i+M ] 2 


(5.8) 


and approaching the origin the radial velocity ^ becomes constant. Therefore, the total affine parameter of the null 
geodesic 


^tot - 


VC r o r°°e 

E , 


(y-To) 2 


dv 


(5.9) 


is found to be finite. Therefore, the spacetime is geodesically incomplete for the null geodesics. 

For the special case c = 0, we have seen that there are two solutions which extend from the origin r = 0 to a finite 
radius and both solutions have divergent curvature invariants at the origin. The first refers to the upper branch of 
e = 1 and the origin corresponds to u —> +oo, while the second refers to the lower branch of e = —1, 0 < a < 3 and 
r = 0 is achieved for u = 0. For the first solution a timelike geodesic with E > \fd can be seen that reaches 
the origin in finite proper time. For the second solution any timelike geodesic reaches the origin in finite proper 
time. Thus, since the solutions are geodesically incomplete and cannot be extended at the divergent origin, they are 
singular. Note however that for the first solution the radial velocity ^ tends to a constant as the particle approaches 
the origin, therefore the radial acceleration is zero. 


VI. DISCUSSION AND CONCLUSIONS 

General branches of new spherically symmetric solutions have been obtained in the context of Quantum Einstein 
Gravity. The source is a “cosmological constant” field that scales at high energies as the asymptotic safety scenario 
suggests, i.e. A (k) oc k 2 . This scaling is also the unique one which is consistent with dimensional analysis without the 
introduction of a new energy scale. The solutions solve the uniquely defined modified Einstein equations which arise 
by adding appropriate covariant kinetic terms to the action in order to ensure the conservation of energy through the 
Bianchi identities. The importance of the presented solutions lies on the fact that they describe the interior small 
distance behaviour of a quantum corrected black hole according to asymptotically safe gravity. It is expected that 
the vacuum quantum modified spherically symmetric solutions in the UV region is dominated by the A field. 

The study of the solutions revealed important new features for the spacetime structure close to the centre of 
spherical symmetry. Since the signature of the metric close to these centres is not known, we have considered both 
cases, the one with the “correct” signs for the temporal-radial metric components (—I- ++) and also the “opposite” 
one with signature (-1-HI-). The “opposite” signature cannot be excluded since the interior Schwarzschild black 


2 Note that the upper branches of cases (ii), (iv) with e = —1, which also extend from the origin to a finite radius and have divergent 
curvature invariants at the origin, are singular since any timelike geodesic can be seen that it needs a finite proper time to reach the 
origin. 
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hole in the same spherical coordinates has also this signature. The physically significant solutions, i.e. those which 
improve or resolve the singularity issues, are separated into two categories. 

The first category consists of one solution with the “correct” signature which starts from the origin r = 0 and is 
accompanied by a curvature singularity at the origin. However, no timelike geodesic can approach the origin since 
any infalling massive particle is reflected at some finite distance. This is due to the repulsive force caused by the 
cosmological constant field along with the other kinetic terms, and this feature is not met by a standard cosmological 
constant. Although the space is timelike geodesically complete (t-completeness), still the massless particles reach the 
origin in finite affine parameter. 

The second category contains various solutions of both signatures which range from a finite radius outwards and 
cannot be extended to involve the centre of spherical symmetry. These solutions have finite all their curvature 
invariants at this minimum radius, so they can also be considered as non-singular physically accepted geometries. 

The presented general solutions which describe the interior of a quantum corrected black hole are valid up to the 
length where the G field starts to become significant. Thus, at that length one should join the presented solutions 
with the quantum corrected one that includes both fields A and G, which consequently for even larger length scales 
should turn into the the classical Schwarzschild black hole. The present work could also be useful in other directions 
that utilize quantum corrected black holes to explain baryogenesis [63], stability issues [fij], or dark energy [65]. 
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